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We study nanoelectromechanical resonators fabricated from suspended flakes of NbSe2 
(thickness'^SO-SO nm) to probe charge density wave (CDW) physics at nanoscale. Variation of 
elastic and electronic properties accompanying the CDW phase transition (around 30 K) are in- 
vestigated simultaneously using the devices as self-sensing heterodyne mixers. Elastic modulus is 
observed to change by 10 per cent, an amount significantly larger than what had been reported ear- 
lier in the case of bulk crystals. We also study the modulation of conductance by electric field effect, 
and examine its relation to the order parameter and the CDW energy gap at the Fermi surface. 



Charge density waves (CDW) [1] result from the cou- 
pling between electrons and phonons, which distorts the 
electron distribution inside the crystal and the lattice 
also undergoes a periodic deformation. The experimen- 
tal techniques used so far to study charge density waves 
include neutron scattering [2], electrical transport mea- 
surements [3] , vibrating reed technique ML angle-resolved 
photoemission spectroscopy (ARPES) [5] amongst oth- 
ers. The CDW modifies both the elastic and electri- 
cal transport properties of the system. Many questions 
regarding the microscopic origin of CDW in transition 
metal dichalcogenides (TMD), which are quasi-2D mate- 
rials, have still not been resolved. In the last few years, 
new insights have emerged about the formation of CDW 
in TMD (issues like Dirac fermion excitations [ij and 
feasibility of Peierls instability 0). In this letter, we 
will focus on the CDW transition in 2H-NbSe2 (a TMD), 
which has elicited a lot of interest for its unique proper- 
ties and is known to support an incommensurate CDW 
at low temperatures (below 35 K). 

With technological advances in small scale device fab- 
rication, CDW in mesoscale and nanoscale systems jll, ^ 
have become the subject of experimental and theoreti- 
cal studies. The isolation of atomically thin NbSe2 lay- 
ers and transport measurements [10] on them have been 
reported. Our study has been conducted on nanoelec- 
tromechanical resonators made from suspended flakes of 
NbSe2 (thicknesS'^30-50 nm, i.e., 45-80 monolayers) in 
doubly clamped geometry. Resonators based on nano- 
electromechanical systems (NEMS) have been realized 
earlier in carbon nanotubes [n| , nanowires 0, 0"3] and 
graphene [14], to name a few materials. They have been 
used in mass sensing, studying the spin-torque effect and 
Casimir force etc. (Ref. [13] and others cited therein) 
Using them to investigate the CDW transition in NbSe2, 
we demonstrate here that such resonators can also act 
as highly sensitive probes for concomitant structural and 
electronic phase transitions at nanoscale. 

To fabricate our devices, we deposit NbSe2 flakes on 
a silicon wafer (coated with 300 nm of insulating Si02) 




FIG. 1: (a) (Left) Resistance (R) as a function of tempera- 
ture (T) for an on-substrate NbSe2 device. There is a feeble 
and broad peak around 35 K (marked by the dashed line), 
characteristic of the CDW transition. (Right) j^^iH^l plot 
showing the CDW feature prominently, (b) Scanning electron 
microscope (SEM) image of a suspended NbSe2 flake. Scale 
bar: 2/im. (c) SEM image showing side view of a different 
suspended flake. Scale bar: l/xm. 



by mechanical exfoliation [lO^. Standard electron beam 
lithography techniques, followed by sputtering Cr and 
Au, are used to establish metallic contacts to these flakes. 
Before metallization of electrodes, the sample is cleaned 
in Ar plasma to remove oxide and organic residue. DC 
electrical transport measurements on such devices show 
a feeble and broad peak in resistance around 35 K due 
to CDW transition (Fig. 1(a)). They undergo supercon- 
ducting transition at 7.2 K. These values are similar to 
observations in bulk crystals. 

The contacted NbSe2 flakes are suspended by etching 
away the Si02 underneath in a dilute buffered HF so- 
lution. The typical trench depth is 180 nm. Fig. 1(b) 
and 1(c) show scanning electron microscope (SEM) im- 
ages of our suspended devices (typical length ~ 2-3 /im). 
To actuate vibrational motion in NbSe2 devices, they are 
used as heterodyne mixers (Fig. 2(a)). The backplane 
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of the sihcon wafer serves as the^gate electrode. An RF 
voltage (amplitude denoted by Vg^) at frequency / and 
a DC bias {Vg^) are simultaneously applied to the gate 

with a bias-tee {Vg^ ^ ^g^)- Due to its capacitive cou- 
pHng to the gate, the NbSe2 flake experiences a vertical 
driving force at the RF frequency / [15]. The source- 
drain voltage (amplitude Vsd) is applied at a frequency 
f^Af^ shifted by A/ with respect to the RF gate sig- 
nal. A/ (~6-20 kHz) is kept constant throughout the 
measurement. The current through the device at the dif- 
ference frequency A/ is called the mixing current (Imix) 



(1) 



where G is the conductance of the NbSe2 flake, q is the 
charge induced by V^^ and ^/ is the amplitude of me- 
chanical vibration (at applied frequency /) along the z- 
direction perpendicular to the plane of the suspended 
flake (see Fig. 2(a)). A and B are factors that depend 
upon applied voltages and the device geometry [18]. 

Since the measurable mixing current depends upon 
how well the conductance is modulated by inducing 
charges (the factor ^ in Eqn. 1), this technique of 
heterodyne mixing had so far been applied to materials 
which show a considerable gating effect (change in con- 
ductance as a function of gate voltage), like carbon nan- 
otubes and semiconductor nanowires. Although NbSe2 
flakes having thickness of tens of nanometers are metal- 
lic, this method works when the resonator is excited in 
the linear regime to such a degree that the resonance 
contribution to Imix is distinguishable against the noise 
floor [l9[. Fig. 2(b) shows a plot of mixing current versus 
frequency for Device 1 at 10 K temperature. The peak 
at 23.9 MHz reflects the mechanical resonance. Quality 
factor of resonance, on fitting to Eqn. 1 (see also Ref. 
^18i]), is 215. In practice, the A/ signal measured by us 
is Imix (given in Eqn. 1) plus a background (of the order 
of hundred pA) arising out of other sources in the circuit. 

At room temperature, when the stress in the flake 
is negligible [20], the natural frequency is given by 

/o=3.56^y^^-^. Here L is the length of the suspended 

flake, E the elastic modulus, / the inertia moment 
(= ^ (width) X (thickness)^), p the mass density (=6.467 
g/cm^) [111 and S the cross-sectional area. From room 
temperature resonant frequency measurements, E of 
most of these NbSe2 flakes of approximately rectangular 
shape are calculated to lie in the range 1.0 — 1.6 x 10^^ 
Pa. On cooHng them down, the metallic electrodes and 
NbSe2 contract, enhancing the tension F and strain e in 
the resonator. In the high tension regime [20j . 
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FIG. 2: (a) Schematic of actuation technique, (b) Measure- 
ment at 10 K, showing mixing current (/mice) as a function of 
driving frequency / for Device 1. (y/^=32 V, A/=12 kHz), 
(c) Colourscale plot of Imix as a function of both frequency 
and DC gate voltage Vg^ for Device 1 at 10 K. The points 
corresponding to mechanical resonance of the device trace out 
a parabolic curve. (A/=12 kHz). 



In Device 1 (length 3.3/im, width 1.7/im, thickness 50 
nm), /o changed from 18.6 MHz at room temperature to 
25.2 MHz at 10 K. Its elastic modulus at room temper- 
ature is 1.0 X 10^^ Pa. Assuming E remains the same 
even at low temperatures and using Eqn. 2, the stress 
(= T/S) at 10 K is calculated to be 1.1 x 10^ Pa (- 1 
kBar) and the strain in the device is 1.1 x 10~^. 

Fig. 2(c) shows a colourscale plot of mixing current as 
a function of frequency and Vg^ (at 10 K). The points 
where mixing current amplitude changes sharply (cor- 
responding to mechanical resonance of the device) are 
distinguishable against the background as the parabola- 
shaped contour. The resonant frequency decreases away 
from V (on both sides) as a function of V^'^. This 
behavior [22| is typical of electrostatically actuated res- 
onators (seen earlier in nanowire [13], nanotube [16] and 
graphene [M] based NEMS devices). 

Temperature scans were conducted by slowly heating 
the devices inside a helium flow cryostat (controlled tem- 
perature sweep at typical rate ~ 75 mK/min in the pres- 
ence of exchange gas to ensure proper thermalization) . 
Fig. 3(a) shows data from Device 2 (thickness 35 nm). 
The striking feature is the huge change in resonant fre- 
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quency around the CDW transition. In the temperature 
range between 28 K and 40 K, change in strain due to 
thermal expansion of Cr/Au electrodes being negligible, 
the resonant frequency can be related to the elastic mod- 
ulus as /o (X \fE (see Eqn. 2). From the observed reso- 
nant frequencies at these two temperatures, the fractional 
change in E across the CDW transition is estimated as 
10 percent |23|]. This signifies a huge degree of softening 
of acoustic phonons in the CDW state, exceedin g by far 
observations in earlier studies on bulk samples 0, US Ell- 
in their vibrating reed study [^], Barmatz et al. reported 
a 0.2 per cent change in elastic modulus and unlike 
their results we do not see a negative peak in £^ as a func- 
tion of T (temperature). The normal - incommensurate 
CDW phase transition in NbSe2 is accepted to be sec- 
ond order. Landau theory calculations suggest a weak 
first order transition in the presence of impurities [26| . 
Lattice distortions and defects can play an extremely vi- 
tal role in such small samples. Impurities can strongly 
pin the incommensurate CDW, modifying its coupling to 
the lattice and this may be one of the reasons for the 
unexpected behavior of elastic modulus observed by us. 

Next we focus on the aspects of the above-mentioned 
data which highlight electronic properties of the CDW. 
The measured resonance feature is prominent in the 
CDW state (below 30 K) but fades out at higher tem- 
peratures (Fig. 3(a)). This rapid change is shown by the 
two plots (at 28 K and 40 K) in Fig. 3(b). The sharp 
peak in Imix corresponding to mechanical resonance al- 
most disappears just after the CDW transition. This 
occurs because ^ reduces drastically while ramping up 
the temperature at the CDW transition. The height of 
the peak in Imix with respect to the background is a 
measure of Its dependence upon the CDW order 
parameter (energy gap at Fermi surface) is borne out 
clearly in Fig. 3(c). Power law fit to ^ - (T^ - T)^ 
gives Tc = 30.6 K and /3 = 0.47±0.1. This estimate of 
P is surprisingly close to 0.5, the critical exponent of the 
order parameter in Ginzburg-Landau theory. When the 
CDW is pinned and gate voltage Vg^ is applied, induced 
charges will modify the carrier density at the edges of the 
'valence' and 'conduction' bands near the Fermi energy, 
leading to a change in the conductance G. (We draw an 
analogy here to the band ga p and electrostatic gating ef- 
fect in a semiconductor [27].) However, the energy gap 
disappears upon transition to the normal state. Con- 
duction becomes metallic, sensitive only to the electron 
concentration at the Fermi energy. In a metal, induced 
charges can not alter the Fermi level significantly because 
of the large density of states. Therefore, ^ (and also 
the Imix signal of resonance) is larger in the presence of 
a CDW than in the normal state. Recently, the existence 
of a CDW energy gap at the Fermi surface of NbSe2 has 
been confirmed by Borisenko et al. using ARPES [Sj. 
The gap at 19.5 K is 2.4 meV, greater than kT (--1.7 
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FIG. 3: (a) Colourscale plot of mixing current, Imix, as a func- 
tion of driving frequency /, scanned over a range of tempera- 
tures for Device 2. The resonant frequency starts to increase 
sharply after 30 K, a signature of transition from the charge 
density wave to normal state. An abrupt change (dashed cir- 
cle) is also noticeable at 18 K. (F/^ = 42 V, A/ = 6 kHz), 
(b) Imix VS. /, at two different temperatures (marked by 
dashed lines in Fig. 3(a), below and above the CDW transi- 
tion temperature, (c) Variation of ^ (arbitrary units) with 
temperature, (d) Finer scan of Imix vs. / and temperature 
around 18 K for Device 1. (F/^ = 28 V, A/ = 12 kHz), (e) 
Imix as a function of / and temperature for Device 1 around 
the CDW transition. (V/^ = - 32 V, A/ = 12 kHz), (f) 
Resonance measurement linescan at 5.4 K (in Device 2) when 
NbSes is a superconductor. (F/^ = 42 V, A/ = 6 kHz). 
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meV). So, the analogy to a semiconductor is justified. 

CDW affects only a small fraction of the Fermi surface 
in NbSe2. There are enough free electrons, even below 
CDW Tc, to screen out any electric field in the bulk. The 
modulation in conductance due to induced charges, giv- 
ing rise to a measurable Imix^ takes place in a few layers 
at the surface of the NbSe2 flake which faces the gate 
electrode. It is interesting to note that we can detect the 
resonance feature prominently even when NbSe2 is su- 
perconducting (Fig. 3(f)). ^ does not show any change 
at or near 7 K, suggesting that the CDW part of the 
Fermi surface is unaffected by superconductivity (a fact 
corroborated by ARPES data [5]) and the CDW energy 
gap exists even in the superconducting state (at least at 
the surface). 

At 18 K, an abrupt decrease in resonant frequency is 
observed for both Devices 1 and 2 (Fig. 3(d) and 3(a)). 
In Fig. 3(d) (data from Device 1), the sharp drop of 280 
kHz corresponds to a 2.1 per cent reduction in E. (In 
Device 2, the change in E at 18 K is 0.7 per cent.) To 
the best of our knowledge, no anomaly has been reported 
earlier in the elastic modulus or thermal expansion coef- 
ficient of bulk NbSe2 crystals around this temperature. 
Apart from lattice imperfections and impurities, the large 
stress (~1 kBar) in our samples can possibly be an im- 
portant factor in determining the elastic properties. Fol- 
lowing the article on structural phase transitions by Re- 
hwald (is), the behavior of E close to 18 K suggests a 
second order phase transition, where the strain couples 
quadratically to the order parameter in the Landau free 
energy. Further experiments are essential. 

Fig. 3(e) shows the variation of E near CDW transi- 
tion of Device 1. However, on a few experimental runs on 
this device, we have observed larger divergent behaviour. 
Occasionally, we have also seen the absence of prominent 
CDW features. These may be signatures of metastability 
or near-critical-point behaviour. This is not well under- 
stood and more detailed studies are necessary to probe 
this phenomenon. 

In conclusion, we have demonstrated a new way to 
look into CDW physics at nanoscale which, as our re- 
sults show, can be different and highly unexpected even 
for well-known materials. This technique can query both 
elastic and electronic properties, and has immense po- 
tential for a wider range of applications to strongly cor- 
related phenomena, including systems undergoing struc- 
tural phase transitions and non-equilibrium properties of 
few vortex dynamics in nanoscale superconductors. 
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